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Abstract 

Despite widespread discussion, the role of van der Waals dispersion forces in wetting remains 
unclear. Here we show that non-local correlations contribute substantially to the water-metal 
bond and that this is an important factor in governing the relative stabilities of wetting layers and 
3D bulk ice. Due to the greater polarizability of the substrate metal atoms, non-local correlations 
between water and the metal exceed those between water molecules within ice. This sheds light on a 
long-standing problem, wherein common density functional theory exchange-correlation functionals 
incorrectly predict that none of the low temperature experimentally characterized ice-like wetting 
layers are thermodynamically stable. 

PACS numbers: 68.43.Bc, 68.43.Fg, 82.65. +r, 
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Water covers almost all solid surfaces under ambient conditions. As such, interfacial 
water is of crucial importance to an endless list of problems in the physical and chemical 
sciences. Examples include heterogenous ice nucleation on aerosol particles of relevance 
to the atmospheric sciences, the flow of confined water through membranes and pores in 
connection with geology and waste water treatment, or the response and reactivity of water 



to electrochemical fields 



A prerequisite to understand these varied phenomena is the 
seemingly simple task of establishing how the water molecules are arranged at the interface, 
i.e., what the water overlayer structure is. Characterizing water overlayer structures is, 
however, a challenging task and despite thousands of publications on the chemical physics 
of water at interfaces only a handful of determinations have been accomplished to date 

n n 

|2|]. These have all been on well-defined atomically smooth substrates - mostly metal 
surfaces - under ultra-high vacuum (UHV) conditions, to which the whole arsenal of surface 
science experimental probes can be applied 

HQ. 

Instrumental in recent determinations 
has been density functional theory (DFT), and indeed there is now an almost symbiotic 
relationship between DFT and scanning tunnelling microscopy (STM), with DFT being 
used to provide structural models with which the experiments can be interpreted. Notable 
examples where DFT has been crucial in unravelling water overlayer structures include water 
in the submonolayer regime on Pd(lll) j^, Cu(llO) [4], and Pt(lll) |5|]. 

Despite the undeniable value of DFT in helping to understand water on metals, there is 
an important and widely discussed problem I3l42l| . Specifically, the structures identified 
with widely used generalized gradient approximation (GGA) functionals such as PBE or 
PW91 are consistently found to be less stable than ice Ih (the most common form of ice) 



191]. This suggests that the intact water overlayers should simply not form under equilibrium 
conditions. This is an uneasy position for the role of DFT in water adsorption studies and it 
casts doubt on the structural characterizations themselves. Since the typical GGAs used in 
such studies do not account properly for van der Waals (vdW) dispersion forces, it would of 
course be interesting and important to know what role vdW forces play in water adsorption. 



Apart from two limited attempts to address this issue 22, l23|], and widespread discussion 



.I3l-l2l|. no detailed study on the role of dispersion forces in realistic water adlayer structures 
on metals has been reported. Thanks to the development of xc functionals which explicitly 
account for dispersion, such as the non-local vdW density functional (vdW-DF) of Dion et 



al. 



2J] and its offspring |25|, it is now possible to tackle this question head on. 



Here we report a study in which the role of dispersion in water metal bonding is examined 
in detail. To this end we apply a newly developed non-local functional to two of the most 
widely studied water-ice adsorption systems, namely water on Cu(llO) and Ru(OOOl). We 
find that non-local correlations contribute substantially to the water-metal bond and that 
this is an important factor in governing the relative stabilities of wetting layers and 3D bulk 
ice. The fact that dispersion forces favor ID and 2D wetting layer structures over 3D bulk 
ice is simply due to the much larger polarizability of the substrate metal atoms. This study 
highlights the key role dispersion plays in wetting and how non-local functionals ameliorate 
the problem common GGAs have with experimentally characterized wetting layers on metal 
surfaces. Since quantitative first principles predictions of wetting and ice growth on metals 
are still beyond reach, some of the outstanding problems are also briefly discussed. 

DFT calculations were performed with a modified version of the VASP 5.2 code 0, [2^, 
which includes our own self-consistent implementation of the non-local van der Waals density 
functionals [28]. Results are reported for PBE [29| and a modified version of vdW-DF, 
referred to as "optB88-vdW" 25|. The difference between the original vdw-DF of Dion et 



al. 



24j | and optB88-vdW is merely in the exchange functional, with the optB88 exchange 



unctional yielding more accurate interaction energies than the original choice of revPBE 



25| . The Cu(llO) and Ru(OOOl) surfaces were modelled with slabs of between 4 and 6 layers 
thickness separated by 14 A of vacuum. Bulk ice was modelled using Hamann's twelve water 
molecule ice cell ^|. Core electrons were treated with the projector augmented wave (PAW) 
method 3l|, whilst the valence electrons {?>d}^As^ and AdJhs^ for Cu and Ru, respectively) 
were expanded in a plane wave basis with a 600 eV cut-off energy. Adsorption energies 
per water molecule were computed from Eads = (-E*°*[H20/M] — _E'*°*[M] — r2i?*°'[H20])/n; 
where i?*°*[H20/M] is the total energy of the n water molecule adsorption system, i?*°*[M] 
is the total energy of the relaxed bare metal slab and i?*°*[H20] is the total energy of an 
isolated gas phase H2O molecule. Likewise the lattice energy of bulk ice is obtained by 
subtracting the total energy of the nH20 ice supercell from the total energy of n gas phase 
H2O molecules. Zero point corrected versions of these quantities were also computed from 
harmonic vibrational frequencies. Some of the energy differences between the structures 
considered here are very small. However, we have carried out extensive convergence tests 
(some of which are reported in Tables SI and SII of the supporting information) which give 
us confidence that our chosen computational set-up is sufficiently close to convergence that 



none of the conclusions reached here would be altered if yet more accurate settings were 
used. 

Let us first consider the wetting of Cu(llO). Below 140 K and for a range of temperatures 
and coverages in the submonolayer regime, water forms ID chains which run perpendicular 
to the Cu(llO) ridges. These chains have been observed by a number of groups {4, 32, 33| 
and although caution must in general be exercised in considering how kinetics dictates water 
film morphology in low temperature UHV experiments, all existing experimental evidence 
suggests that the ID chains are the thermodjTiamically stable phase for water on Cu(llO) at 
low coverage. The structure of the chains has recently been characterized and shown to be 
a periodic arrangement of water pentagons on the basis of STM, infra-red spectroscopy, and 
DFT [4]. Hence the ID chains provide an excellent opportunity to understand the relative 
stability of bulk ice with a well-defined wetting layer and to explore the role of dispersion 
forces in wetting on metals. 

The structure of the water chains is shown in Fig. [D^a). They are comprised of a 
face sharing arrangement of water pentagons with two types of water molecules: fiat-lying 
molecules bonded directly to Cu surface atoms along the ridges and upright molecules over 
the troughs in the (110) surface which interact relatively weakly with the substrate. Fig. [2] 
summarizes the adsorption energies obtained, along with the corresponding lattice energies 
of bulk ice. Results with and without zero point energies are reported. The first observation 
is that PBE predicts bulk ice to be more stable than the pentagon-based wetting layer. 
Specifically, the energy of the adsorbed structure is 92 meV/H20 less stable than the lattice 
energy of bulk ice. Accounting for the zero point energy difference between the overlayer 
and ice is not sufficient to correct the ordering of the two states, although it does push things 
in the right direction by reducing the preference for bulk ice to 63 meV/H20. This incorrect 
preference for ice is precisely the problem that has been observed many times before for 
pure water overlayers computed with PBE and PW91 . Moving to the optB88-vdW 

functional, we find that it does indeed give the correct energy ordering, predicting that the 
overlayer is marginally (13 meV/H20) more stable than bulk ice at the total energy level 
and 39 meV/H20 more stable when zero point energies are accounted for. Thus a stable 
wetting layer is found in agreement with the experimental observations. The structure of 
the overlayers obtained with PBE and optB88-vdW are very similar (all hydrogen bonds 
between water molecules and water-metal distances are within 0.02 A), which indicates 




FIG. 1. Top views of water overlayers: (a) a pentagonal chain on Cu(llO) [4]; (b) H-down bilayer 

I ^ 

on Cu(llO); (c) extended chains on Ru(OOOl) with flat-lying and H-down molecules [15]; and (d) 
H-down bilayer on Ru(OOOl). Black, red, dark gray, and light gray circles represent H, O, and Cu 
(Ru) in the first and second layer, respectively. For molecules in the H-down configuration the H 
atoms beneath the oxygens are not visible. The unit cells used are indicated. 

that whilst dispersion is important for the absolute adsorption energy of water it is of less 
importance to the atomic structure. 

As a consistency check, we compared the adsorption energy of the pentagonal overlayer 
with extended 2D bilayers (both the "H-down" , Fig. WO^), and "H-up" bilayers considered for 
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FIG. 2. Adsorption energies of water overlayers on Cu(llO) and Ru(OOOl) compared to the lattice 
energy of bulk ice Ih using PBE and optB88-vdW. The adsorption energy is also decomposed in to 
water-water and water-metal bonding. The thin centered black bars and associated numbers are 
the adsorption/lattice energies corrected for zero point energies. 

this system before [l^) and alternative ID chain models Irrespective of the functional 
used the pentagonal chain structure remains the most stable overlayer. Therefore, the 
relative energies of the various overlayers are not strongly affected by dispersion forces, 
which explains why PBE and PW91 have had success in predicting water structures on 
metals. Finally, explicit consideration of finite temperature vibrational free energy effects 
does not alter the conclusions drawn here and when these effects are accounted for the 
pentagons remain more stable than bulk ice (Fig. SII). 

Regarding now the wetting of Ru(OOOl), Hodgson and co-workers have shown that be- 
tween 140 and 160 K a 2D wetting layer of intact water molecules forms, which upon heat- 
mg or exposure ,Q. At 0.67 



TABLE I. Different contributions to the optB88-vdW adsorption energy (-Eads) of the pentagons on 
Cu(llO) and the extended chain structure on Ru(OOOl): water-water bonding {E^^^~^'^'~'), water- 
metal bonding {E^^^~^), non-local correlation contribution to the adsorption energy {E'^^^), non- 
local correlation contribution to water- water bonding (£'^1^^^^'-''°''^), and non-local correlation 
contribution to water- metal bonding iE^^^~^'^^'^). Zero point energy corrected values are given 
in parenthesis. All values in meV/H20. 
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monolayers the intact water overlayer exhibits a sharp \/3 diffraction pattern in low-energy 
electron diffraction (LEED) but a very low specular reflectivity in helium atom scattering 
(HAS). Although the precise structure of the overlayer has not been resolved, a compelling 
model that explains both the order observed in LEED and the apparent disorder from HAS 
has been put forward. The overlayer involves water adsorbed near Ru atop sites in a hy- 
drogen bonded honeycomb network, containing chains of fiat-lying and chains of H-down 
bonded water in a hexagonal superstructure. Here, we have considered this extended chains 
model in a 2a/3 x ^/3 unit cell (Fig. lH^c)), as well as intact H-down (Fig. [H^d)) and H-up 
a/3 X a/3 bilayers. With PBE the chains are more stable than the bilayers but are less 



stable than bulk ice (Fig. 2 and Table SH) 3J]. Therefore, as with water on Cu(llO), with 



PBE one would not expect the formation of an extended overlayer. However, when vdW 
interactions are accounted for the situation is reversed. In particular, optB88-vdW predicts 
that the extended chain structure is indeed 72 meV/H20 (97 meV/H20 if zero point en- 
ergy is included) more stable than bulk ice. These results are in apparent agreement with 
experiment, where wetting is observed for an intact water overlayer. Again the structures 
obtained from PBE and optB88-vdW do not differ to any great extent (all distances are 
within 0.02 A). 

We now discuss why dispersion improves the relative energies of the overlayers and bulk 
ice. A decomposition of the total adsorption energy in to water-water and water-metal 
contributions proves useful (see the supporting information for details). With PBE the 
total adsorption energy is split equally between water-water and water-metal bonding (Fig. 



[2]). When we switch to optB88-vdW, the adsorption energy increases by about 150 and 
180 meV for Cu and Ru, respectively. Crucially this increase comes almost exclusively 
from water- metal bonding, with essentially no change in water- water bonding (Fig. 2). 
Further analysis shows that this increase in water-metal bonding is largely due to non-local 
correlation. Indeed, on both Cu and Ru the non-local contribution to the water-metal 
bonding exceeds that for the water- water bonding by a factor of 2 to 3 (Table [T]). It can 
also be seen from Table [T] that the non-local contribution to the water- metal bonding on Ru 
is larger than that on Cu, which explains why dispersion has a larger impact on Ru than 
on Cu. Although, it is somewhat counterintuitive that dispersion favours the formation of 
ID and 2D structures over 3D bulk ice, our analysis is consistent with the fact that Ru 
and Cu atoms have much larger polarizabilities than O and H. Specifically, the relative Cq 
coefficients from time-dependent DFT calculations 35| of Ru : Cu : O : H are 79 : 44 : 2 : 
1, respectively. 

Before concluding, we stress that although progress has been made we are still some 
way from quantitative first principles predictions of wetting and related phenomena such 
as heterogeneous ice nucleation rates. In particular, here we have obtained an improved 
description of the relative energies of water-ice overlayers and bulk ice, but this is achieved at 
the expense of a slightly worse absolute lattice energy for ice. The experimental lattice energy 
for bulk ice Ih is 610 meV/H20 16|, our PEE value is -663 meV/H20, and our optB88- 
vdW value is -711 meV/H20. PBE outperforms optB88-vdW on the lattice energy because 
the missing dispersion is compensated for by the tendency of PBE to slightly overestimate 



H-bonds between water molecules which are almost linear |36|. optB88-vdW, on the other 
hand, is known to overestimate the vdW correlation as the number of H-bonds between 
water molecules increases 



25|. Associated with this overbinding is a lattice constant that is 
too small: the experimental volume is 32.1 A'^/H20, the PBE value is 30.4 A^/Il20, and the 
optB88-vdW value is 29.8 A^/H20. The issue of the lattice constant of bulk ice is particularly 
relevant to heterogeneous ice nucleation, where epitaxial match between the substrate and 
ice has long been considered a key factor in the nucleating ability of a substrate. It remains 
to be seen if non-local vdW functionals in their current form are capable of simultaneously 
addressing all of these issues. Analysis of the underlying errors in the exchange and (non- 



local) correlation components might help to shed light on this 37|]. Although, in the longer 



term the application of the random phase approximation within the adiabatic connection 
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fluctuation dissipation theorem or quantum Monte Carlo to ice-like overlayers on metals 
would be welcomed. 

In conclusion, we have considered the role of dispersion in water-metal bonding. Analysis 
of the relative contributions of water- water and water-metal bonding, to two wetting layers 
of contemporary importance, highlights the critical role dispersion plays in the wetting of 
metals by water. Given that DFT plays a central role in interpreting and understanding 
the most well-defined experimental studies of wetting, it is satisfying that progress with the 
long-standing wetting problem has been made. Although dispersion is of minor importance 
for many questions of structure and bonding, we have shown here that the dispersion forces 
between water and metals are sufficiently large to favor formation of ID and 2D wetting 
structures over 3D bulk ice. This arises from the much larger polarizability of the metal 
atoms compared to oxygen and hydrogen, and so it is likely that dispersion will be of 
importance to water on metals in general and not just the exemplar systems considered 
here. 
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